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One of the major goals of ultrarelativistic nuclear collision studies is to obtain informa- 
tion on the QCD phase diagram [lj. A promising approach is the investigation of hadron 
production. Hadron yields measured in central heavy ion collisions from AGS up to RHIC 
energies can be described very well (see [5] and refs. therein) within a hadro-chemical 
equilibrium model. In our approach the only parameters are the chemical freeze-out tem- 
perature T, the baryo-chemical potential //& and the fireball volume V (for a review see 
0). 

The main result of these investigations was that the extracted temperature values rise 
rather sharply from low energies on towards ^s^n —10 GeV and reach afterwards con- 
stant values near T=160 MeV, while the baryochemical potential decreases smoothly as 
a function of energy. The limiting temperature [3] behavior suggests a connection to the 
phase boundary and it was, indeed, argued [4] that the quark-hadron phase transition 
drives the equilibration dynamically, at least for SPS energies and above. Considering 
also the results obtained for elementary collisions, where similar analyses of hadron mul- 
tiplicities, albeit with several additional, non-statistical parameters (see [5116] and refs. 
therein), yield also temperature values in the range of 160 MeV, alternative interpreta- 
tions were put forward. These include conjectures that the thermodynamical state is not 
reached by dynamical equilibration among constituents but rather is a generic fingerprint 
of hadronization [7118] . or is a feature of the excited QCD vacuum [9]. 

While in general all hadron yields are described rather quantitatively, a notable excep- 
tion was up-to-now the energy dependence of the K + /tc + ratio which exhibits a rather 
marked maximum, "the horn" [TO], near ^s^n — 10 GeV [TTJ. Predicted first within a 
model of quark-gluon plasma (QGP) formation [TO], the existence of such a maximum 
was also predicted [12] within the framework of the statistical model, but the observed 
rather sharp structure could not be reproduced [13] (see also the discussion in [TT]). As a 
consequence, the horn structure is taken in [11] as experimental evidence for the onset of 
deconfinement and QGP formation, and as support for the predictions of [TO]. We have 
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recently shown [TJ] that, employing an improved hadronic mass spectrum, in which the a 
meson and many higher-lying resonances are included, leads to a sharpening the structure 
in the K + /tt + ratio, as will be shown below. 
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Figure 1. The energy dependence of the 
temperature and baryon chemical potential 
at chemical freeze-out. The results of the 
new fits [13] are compared to the values ob- 
tained in our earlier study |13j . The lines are 
parametrizations for T and fib (see text). 



Figure 2. Energy dependence of the relative 
production ratios K + /n + and A/tt~. With 
the dotted line we show for the K + /tc + ratio 
an estimate of the effect of higher mass res- 
onances. The dashed lines show the energy 
dependence of T and fib- 



In Fig. [T] we present the energy dependence of T and fib in comparison to results of 
other similar analyses P^|16II17P^] . In [TS] the strangeness suppression factor 7 S is an 
additional fit parameter; in [T7j an inhomogeneous freeze-out scenario is modeled with 
spreads in T and fib as extra fit parameters; the approach of [19], employing a full set of 
non-equilibrium fit parameters, leads to rather different results compared to those shown 
here. We have parametrized our values of T and function of ^Snn with the 

following expressions: 

1 1 SCR 
T = T lim , fc b [MeV] = (1) 

l + exp(2.60-ln( v /s iV7 v(GeV))/0.45) 1 + 0.28QJ s NN {GeV) 
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with the "limiting" temperature T^ m =164 MeV. 

We employ the above parametrizations of T and fib to investigate the energy depen- 
dence of the production yields of K + and A hadrons relative to pions, shown in Fig. [2J 
The K + /tt + ratio shows a rather pronounced maximum at a beam energy of 30 AGeV 
and the data are well reproduced by the model calculations. In the thermal model 
this maximum occurs naturally at v /saT/v —8 GeV [12]. It is due to the counteracting 
effects of the steep rise and saturation of T and the strong monotonous decrease in fi h - 
The competing effects are most prominently reflected in the energy dependence of the A 
hyperon to pion ratio (lower panel of Fig. [2]), which shows a pronounced maximum at 
t/snn —5 GeV. This is reflected in the K + /tt + ratio somewhat less directly; it appears 
mainly as a consequence of strangeness neutrality, assumed in our calculations. 

The model describes the K + /tt + data very well over the full energy range, as a con- 
sequence of the inclusion in the code of the high-mass resonances and of the a meson, 
while our earlier calculations [13] were overpredicting the SPS data. At RHIC energies, 
the quality of the present fits is essentially unchanged compared to [13], as also the data 
have changed somewhat. The model also describes accurately the h/ix~ measurements as 
well as those for other hyperons |14j . 
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Figure 3. The phase diagram of 
strongly interacting matter. The 
points represent the results of the 
thermal fits. For the SPS beam en- 
ergy of 40 AGeV (fi b ~400 MeV) we 
show both midrapidity (dN/dy) and 
full phase space (47r) fit results. The 
phase boundary and critical point 
from lattice QCD (LQCD) calcu- 
lations [20] is shown together with 
freeze-out curves for a hadron gas 
at constant baryon density (baryons 
and anti-baryons) and energy den- 
sity. The full triangle indicates 
the location of ground state nuclear 
matter (atomic nuclei). 



In Fig. [3] we show the result of our fits in the phase diagram of strongly interacting 
matter (for a recent review see [22] ) . Our results strongly imply that hadronic observables 
near and above the horn structure at a beam energy of 30-40 AGeV {fib —400 MeV), coin- 
ciding with the approach to saturation in T, provide a link to the QCD phase transition. 
Open questions are whether the chemical freeze-out curve for larger values of fib actually 
traces the QCD phase boundary or whether chemical freeze-out in this energy range is 
influenced by exotic new phases such as have been predicted in [21J. 

In summary, our recent results [H] demonstrate that by inclusion of the a meson and 
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many higher mass resonances into the resonance spectrum employed in the statistical 
model calculations an improved description is obtained of hadron production in central 
nucleus-nucleus collisions at ultra-relativistic energies. A dramatic improvement is visible 
for the K + /tt + ratio, which is now well described at all energies. The "horn" finds 
herewith a natural explanation which is, however, deeply rooted in and connected to 
detailed features of the hadronic mass spectrum which leads to a limiting temperature 
and contains the QCD phase transition [3]. 

We acknowledge the support from the Alliance Program of the Helmholtz-Gemeinschaft. 
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